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Abstract 

We consider the one-dimensional XXX spin 1/2 Heisenberg antiferromagnet at zero temperature 
and zero magnetic field. We are interested in a probability of formation of a ferromagnetic string 
P{n) in the antiferromagnetic ground-state. We call it emptiness formation probability [EFP]. 
We suggest a new technique for computation of the EFP in the inhomogeneous case. It is 
based on the quantum Knizhnik-Zamolodchikov equation [qKZ]. We calculate EFP for n < 6 
for inhomogeneous case. The homogeneous limit confirms our hypothesis about the relation of 
quantum correlations and number theory. We also make a conjecture about a structure of EFP 
for arbitrary n. 
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1. Introduction 



The Hamiltonian of the XXX Heisenberg spin chain can be written hke this 



N 

H= i: + ^r<i + o^K+i -1) (1-1) 

i=-N 

Here 2N+ 1 is the length of the lattice and , af, are Pauli matrices. We consider thermodynamic 
limit [ goes to infinity] . The sign in front of the Hamiltonian indicates that we are dealing with 
the antiferromagnetic case. We also imply periodic boundary conditions. The model was solved by 
Bethe in 1931, see [[l|. The ground state was constructed by by Hulthen, see p|. We shall denote 
the ground state in the thermodynamic limit by |GS) . The emptiness formation probability (EFP) 
for the XXX model is defined as follows: 

n 

P(n) = (GS|n^.|GS), (1.2) 

i=i 

here Pj = Sj + ^ is the projector on the state with the spin up in the j-th lattice site. The integer n 
has a meaning of a length of a continuous ferromagnetic string. P{n) is a probability that this string 
can appear in the antiferromagnetic ground-state. The importance of EFP was emphasized in 0. 

P{n) can be represented as n-multiple integral. The integral representation follows from the 
work of RIMS group. RIMS approach is based on vertex operators and bosonic representation of 
infinite-dimensional quantum algebras, see |^. The explicit formula for P{n) in XXX limit was 
obtained in 0: 

= n / ^-(^1' • • • ' Tn{Xl, ...,Xn) (1.3) 

where 

jj , . _ ni<fc<,<nSinh7r(A, -Afc) 

f/„(Ai,...,A„) - TT n-=isinh'^vrA, ^ ^ ^ 

and 

lll<fc<j<nl^j ^ — 'I') 

The contour C goes parallel to the real axis with the imaginary part confined between and —i 
for each integral. 

In papers [|, ^ [1^] we evaluated the integrals for n < 5. We discovered that these P(n) can 
be expressed in terms of the values of Riemann zeta function at odd arguments, log 2 and rational 
coefficients. We conjectured that this is a general property for all P{n). In this paper we proved the 
property for P(6). We think that all correlation functions 

(GS|<< ...<|GS) (1.6) 

also have this property. Asymptotic behavior of P{n) for large n was studied in the papers |l8|, [p^ , 
[|13| and [Q. The technique of calculation of these integrals, described in the paper worked 
for n = 2, 3, 4. In the paper we calculated P{5) by means of this technique. However, these 
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computations are so complicated that it is problematic to generalize them to the case n = Q. So we 
start looking for indirect methods of evaluation of the integrals (1.^. 

It appeared to be useful to consider inhomogeneous case. In this case there are more free param- 
eters. We call them inhomogeneity parameters and denote hj Zi, . . . , Zn- The EFP in the inhomoge- 
neous case we shall denote by P„(zi, . . . , Zn). Let us remind to the reader that inhomogeneous models 
were used for evaluation of correlation functions from the very beginning. For the massive regime 
of the XXZ model the vertex operator approach was developed in @, ^. It allowed to express the 
correlation functions in terms of the trace functions. Special combinations of these trace functions 
satisfy the quantum Knizhnik-Zamolodchikov equation (qKZ), see [11, 1^. Later in paper Miwa 
and Jimbo suggested that the correlation functions in the gap-less regime directly satisfy the qKZ 
equation. Since the XXX model belongs to the gap-less regime we shall use qKZ for evaluation of 
EFP in the inhomogeneous case. We suggest a general ansatz for P„(zi, . . . , 2;„), see (3.20). This 
constitute a new method for computation of the EFP. On the other hand, it is easy to generalize 
the technique explained in ^ to the inhomogeneous case and calculate the EFP directly [for short 
ferromagnetic strings]. When we can compare results, they coincide. 

The paper is organized as follows. In the Section 2 we discuss the relation of EFP to the qKZ 
and derive three important properties of the EFP in the inhomogeneous case. In Section 3 we apply 
a generalization of the technique described in [^ to the inhomogeneous case and compute P„ directly 
for n < 4. Then we check that these P„ satisfy all the properties, which follow from qKZ. This 
helps us to formulate a general ansatz for P„ in the inhomogeneous case. Further, we suggest a 
new way of computing P„. One can use the ansatz and general properties of P„ [which follow from 
the qKZ]. In this way we get the explicit expressions for P5 and Pq in the inhomogeneous case. 
In Section 4 we discuss the homogeneous limit of P{n) for n < 6. In particular, when n < 5 we 
reproduce our previous results, obtained in ||, ^ |T0[. We also get the analytic expression for P(6) 
in the homogeneous limit. Having this answer we can compare it with the numerical value for P(6) 



obtained by the DMRG method in \TU\. We also discuss the structure of EFP in the homogeneous 
case and offer some plausible conjectures. In the last Section 6 we discuss the results and outline 
some possible ways of a further progress. 



2. The EFP in the inhomogeneous case and the qKZ equa- 
tion. 



We believe that consideration of the inhomogeneous case instead of the homogeneous one can give 
us a new information about the EFP and other correlation functions. 

A method of calculation of correlation functions, which we use was found in the papers [^, ^. It 
is based on theory of infinite-dimensional quantized algebras and vertex operators. We shall need 
elements of this method. Let us introduce some notations. 

We use the R-matrix: 



P(A) 



X + ni 











-A 






-X 








X + ni J 



(2.1) 
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where 



Ro{\) 



r (^] r f i - ^ 



V 27ri y V 2 27ri y 

notice that 

Ro{X)Roi-X) = 1 

This R-matrix appears in the rational hmit from XXZ R-matrix, it is related to usual XXX R-matrix 

by 

R{X) = {a'^I)Rxxx{X){I^(T') 

(similar transformation is needed when obtaining form factors of S'f/(2)-invariant Thirring model 
from SG ones |TB[)- 

This R-matrix (2.|l|) satisfies the equation: 



R{-7ii) 



/ 

110 

110 

V y 



(2.2) 



Following [0 we introduce functions gn which satisfy the qKZ equations on level —4 |T^. We 
write the qKZ equations in their original form |]T3|, |T^ which takes into account explicitly symmetry: 

fi'n(Al, ■ ■ ■ , Xj + l, Xj, ■ ■ ■ , X2n)ei,-,t'.^^,e'^,-,e2„ = (2-3) 
= R{Xj - Aj + i) r^ll gn{Xl, ■ ■ ■ , Xj, Xj+i, ■ ■ ■ , X2n)ei,---,ej,ej+i,-,t2n 

dniXi, ■ ■ ■ , X2n-1, X2n + '^'^i) ei,- ,t2u-\,e2n = 9n{X2n-, -^1, " ' ' ) A2n-l)e2n,ei,---,e2n-l (2-4) 

Solutions to these equations are meromorphic functions with possible singularities at the points 

^[\. - Afc) = tt/, le Z\0 

For application to the correlation functions we are interested in a particular solution described 
in details in 0, §]. Detailed study of this solution will be performed in a future publication, in the 
present paper we need only limited information about it. First, the gn is regular at ^{Xj — Xk) = ±7r. 
Moreover, much can be said about its values at these points pi : 



9n{Xi, ■ ■ ■ , Aj_i, Xj, Xj — iri, Xj+2, ■ ■ ■ , A2n)ei,-.-,ej_i,e^,e^+i,ej+2,---,e2n — 
= ^tj-tj + i 9n-l{Xi, ■ ■ ■ , Aj_i, Aj+2, ■ ■ ■ , A2n)£i,-,ej_i,ej+2, -,e2n 

Together with the symmetry (2.3), and (2.|D this equation implies 

X! 9n{Xl, ■ ■ ■ , Aj_i, Xj, Xj + iri, Xj+2, ■ ■ ■ , X2n)ei,-,ej.i,ej,ej+i,ej+2,-,e2„ 
= fl'n-l(Al, ■ ■ ■ , Aj_l, Aj+2, ■ ■ ■ , X2n)ei,-,ej_i,ej+2,-,e2n 

The emptiness formation probability P„ is related to gn as follows 

Pnizi, ■■■ ,Zn) = gn (tTZi, ■ ■ ■ , TTZn, rr{Zn + i) , ■ ■ ■ , Tl{zi + i)) , . , 



(2.5) 



(2.6) 



(2.7) 
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Now we want to establish some general properties of P„ following from (2.3), (2.^), (2.5), (2.6). 

1. Symmetry. The function _P„(zi, ■ ■ ■ , Zn) is symmetric. 
Proof. Obviously, it is enough to show that 

This identity follows from (2.3) and from the fact that the R-matrix acts diagonally on the indices 

-, - and +, +: 

Pni' ' ' ) ^j) ' ' ') 

= 9ni---, TTZj, TTZj+i, ■ ■ ■ , 1T{Zj+i + i), Tc{Zj + i) ■ ■ •)...__,...,+,+,... = 

= Roiii{zj+i - Zj))Ro{n{zj - zj+i)) 

X fi-n (■ ■ ■ , TTZj, ■ ■ ■ , n{Zj + i), Tc{Zj+i + i) ■ ■ ■)...__,...,+,+,... = 

= Pn{' ■ ■ ; Zji ' ' ■) 

QED 

2. Vanishing. The function P„(zi, ■ ■ ■ , -2„) vanishes when z^ = Zj + i. 

Proof. Due to the previous property it is sufficient to consider the case k = 1, j = 2. Let us put, 
first, zi = z + i, Z2 = z', then we shall take the limit z — > z'. 

= gn (vr(z + i),7cz\ ■ ■ ■ , 7rz„, 7i{zn + z), ■ ■ ■ , 7c{z' + i),n{z + 2i)) 



gn {tTZ, 77(2; + i), TTZ', ■ ■ ■ , TTZn, Vr(z„ + «),■■■, T[{z' + i), [ 



'+,-,-,■■■,-,+,■■■,+ 

where we used (2.§). Consider the limit z — > z'. As it has been explained singularities do not occur 
for Q{Xj — Xk) = ±7r. Moreover, the final result contains the fragment 

gni---,TT{z + i),'KZ---)...__^... 

which implies that the result vanishes due to (2.5). Because of absence of singularities this zero does 
not interfere with any pole, so 

P„(z + z,^,---,z„) = (2.8) 

QED 

3. Normalization. The following asymptotic holds for zi — > 00 along the real axis: 

Pn{Zl, Z2, - ■ ■ , Zn) — ^ \Pn-l{z2-, " " " , Z^) (2.9) 

Proof. One more property of the solution gn will be important for us. Using the integral formula 
from one can show that 5'n(Ai, • • • , X2n-2, A, A + ttz) behaves as 0(1) when X ^ oo + in where k 
is a finite number. The leading term of asymptotic does not depend on k. We shall use notation: 

fi'n(Al, ■ ■ ■ , A2n-2, A, A + T^i) e-i,-,t2n^2,e2n^i,e2u ^ ^n(Al, " " " , A2n-2)ei, -,£2,1-2 ; £2n-l,e2n (2.10) 

The function cjn possesses important property of symmetry with respect to last two indices because 

5'n(Ai, ■ ■ ■ , A2„-2, A, A + 7ri)ei, 

= gn{X - ni, Ai, ■ ■ ■ , A2n-2, A), 

= P(Ai - A + TrOfi;^--^ ■ ■ ■ P(A2„-2 - A + my.r^X.T 
X 9n{Xi, ■ ■ ■ , A2„-2, A - TTi, A),'^,...,4^_2,4^,,2n-i 

— *• Sign(e2n-1, ^2n) dn{Xi, ■ ■ ■ , A2n-2)£i,- -,e2n-2 ; e2n,£2n-l (2-11) 



■,e2n-2,e2n-l,e2n 



/^2n^l,---,^2n-2,<^2n-l 
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with the sign function 



signfc,.,) = { / l^Z\ (2.12) 



where we used the asymptotic of the R-matrix 

H) ■ diag(l, -1, -1, 1) 

From the equation (2.11) we conclude that 

dni^l, ■ ■ ■ , ^2n~2)ei,-,t2n-2 ; e2„-i,e2„ =0 if e2n-l = G2n 
dni^l, ■ ■ ■ , A2n-2)ei,-,e2n-2 ; +,- = 9n{^l, " " " , ^2n-2) ei,-,e2„-2 ; -,+ (2.13) 

The relation (2.|l3|) allows to calculate ^n(Ai, ■ ■ ■ , X2n-2)ti,---,e2„-2 ; +,-• 

5'ri-l(Al, ■ ■ ■ , A2n-2)ei, -,£2„-2 = 
= dni^l, ■ ■ ■ , A2n-2)ei,-,e2n-2 ; +,- + ^"("^l; ' ' ' ) X2n-2)ti,-,t2n-2 ; -,+ = 
= 2 ^n(Ai, ■ ■ ■ , A2n-2)ei, -,e2u-2 ; + ,- 

Now the normalization (2.^) follows from 

= Qn {l^Zi, 'KZ2, ■ ■ ■ , TTZn, n{Zn + i) , ■ ' ' , ^{^2 + 0? ^(^1 + 0)- -,...-,+,...,+,+ = 

= gn {nz2, nzn, rr{zn + n{z2 + i), it{zi + i), tt{zi + 2i))_^... _ 

9n (vr^2, ■ ■ ■ , 7rz„, Ti{zn + 0> ■ ■ ■ , 7r(z2 + «))_,...,_,+,...+ . = 

= \ Pn-l{z2,---,Zn) (2.14) 

where the formula (2. PUD was used. QED 
The above three properties will be very useful for further consideration. 

3. Explicit expressions for EFP. 

In paper Jimbo and Miwa suggested an integral representation as a solution to the qKZ (2.3) - (2.|) 
which also satisfies the property (2.5). Using the relation (2.|^) one gets the integral representation 
which is the direct generalization of the formula (1.0) to the inhomogeneous case 

n(n+i) sinh {7i{zk - Zj)) pi/2+oo r-i/2+oo dX^ sinh(7r(Aj - A^)) 

, , , Tiizh — z^) J-i/2-oo Im J-i/2-oo t:^, Ai — At — z 



'^{Zk - Zj) J-i/2~oo 2TXi"' J-i/2-oo 27li j^^. {\j - Afc - z) 

n •=i(m=UA, -zk-t) nL,+i(A, - z,)) 



(3.1) 



n-=inLisinh(7r(A,-Zfc)) 

where for the moment we will consider the inhomogeneity parameters zj as distinct real numbers. 
These parameters also may be considered as arbitrary complex numbers, but the integration contours 
in this case should be taken in such a way that they separate the singularities of the integrand in the 
same manner as the contours in the above formula (3.^) for the real values Zj. 
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For evaluation of the integral (3.|I]) we can use the technique described in 0. We obtain for the 
first four values of P„(zi, . . . , Zn) '■ 

Pi = l (3.2) 



^2 = ^ + ^2) G{zi - Z2) 



(3.3) 



-P3 = ^ + ^3,1(^1, Z2\z-i)G{zi - Z2) + ^3,1(^1, 23^2)^(^1 - z^) + ^3,1(2:2, Z'i\zi)G{z2 - 2:3) (3.4) 

-P4 = 7 + ^4,1(2:1, Z2\z^, Zi)G{zi - Z2) + ^4,1(^1, Z2.\z2, za) G{zi - z^) + v44,i(zi, 2:412:2, 2:3) G{zi - 2:4) + 
5 

M,l{z2, Z'ilZi, Z4) G{Z2 - Z3) + ^4,1(^2, Zi\zi, 2:3) ^(2:2 - 2:4) + ^4,1(2:3, 2:4, \zi, Z2) G{Z3 - 2:4) + 
^2(^1, Z2\Z3, Zi)G{Zi - Z2)G{Z'i - 2:4) + ^4,2(^1, ^3^2, ^4)^(2:1 - 2;3)G(2;2 - 2:4) + 

-44,2(^1, ^4|^2, 2^3)^(2:1 - 2;4)G(2;2 - 2:3) 



where 



G{x) = {x^ + l){p{ix) + Pi-ix)) = 2 ^(-1)^ ■ k ■ 



and the functions 



k=l 



+ 1 



/3(rc) 



i«„i±£)-„£)).ghi£ 



d 

dx 



In r(x) 



(3.5) 
(3.6) 



with the polynomials 



^2,1(^1,^2) = Q2,l{Zl,Z2) 

^3,1(^1. ^2^3) 



^3,1(2:1, 2:212:3) 
^4,2(2:1, 2:212:3, 2:4) 
^A,l{Zl, Z2\Z3, Z4) 



Z13Z23 
Qi,2{Zl-, Z2\Z3, 2:4) 

2;i32;i42;232;24 
^4,1(2:1, 2:21^3, ^4) 

2;i32;i42:232;24 



Q2a{zi-, Z2) 



1 
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^3,1 (2:1, 2:212:3) = —(2:132:23-1) 
Q4 2(2:1, 2:212:3, 2:4) = ^ I (2:13^24 - 1) (2:142:23 - 1) + ^(2:12 + ^)(^34 + ^) + ^ 

(^4 1(^1,2:212:3, 2:4) = 2^4 2(2:1,2:21^3,2:4) - ;^(2:i2 + 4)(z|4 + 1) 



(3.7) 

(3.8) 
(3.9) 
(3.10) 
(3.11) 
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and by definition 

The function G{x) (3.|]) is actually real for real x and even 

G{-x) = G{x) (3.12) 

Besides, for small x one has an expansion 

oo 

G(x) = -2(1 + x^) ^(-1) V=Ca(2A; + 1) (3.13) 

k=0 

with the alternating zeta series defined in (4.|l]). In particular, 

G(0) = -21n2 (3.14) 

The following limit is well defined 



G{±oo) = -- (3.15) 



We shall also need the following relation 

G{x -i) = a{x) + -f{x)G{x) (3.16) 

where 

(x — 2i) , , x(x — 2i) 

^ ' [x-i] ^ ' (x + i){x -i) ^ ' 

In particular, 

G{±%) = -2 (3.18) 

Looking at the answers (3.0-3.^) we can conclude that they satisfy the general three properties 
from the previous section and one more property, namely, the translational invariance: 

Pn{Zi + a, . . . ,Zn + a) = Pn{Zi, . . . ,Zn) (3.19) 

which actually follows from the integral representation (3.|l]). This property means that P„ depends 
only on differences of n parameters Zj. 
One can also suggest a general ansatz : 

[f] / 
Pni,Zi, . . . , Zn ) = ^{And{zi, Z2\z3, Z4\... \z21-1, Z2i\z2i+i . . . Zn)\\ G{z2j-i - Z2j) + permutatlous} 

1=0 j=l 

(3.20) 

where Anj^zi, Z2\z3, Z4\ . . .\z21-1, Z2i\z2i+i ■ ■ ■ Zn) are some rational functions which depend only on 
differences Zj—Zk- We also imply that these functions are symmetric under independent transpositions 
Zi Z2, z^ ^ Z4, . . . , Z21-1 Z21 and any permutation of the residual variables Z21+1 ■ ■ ■ z^. Also we 
require a symmetry of this function under transposition of any pair Z2r-iZ2r ^ Z2s-iZ2s for r,s < I. 
It is implied that the permutations in (3j20|) do not involve those sets of variables for which the 
function An^i is already symmetric. We also take 

AnfliZi,. . . ,Zn) = (3.21) 

n + 1 
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and a 

A r \ \ \ I _ Qn,l{Zi, Z2\Z3, . . .\Z2i-l, Z2l\z2l+l . . . Zn) /o oo^ 

^n,il^l) ^2p3; • • • p2Z-l, Z2l\Z2l+i . . . Zn) — — r [6.AA) 

lll<j<k<n[Zj — Zk) 

where Qn,i are some polynomials of maximum power n — 1 for each variable Zk. Moreover we expect 
that these polynomials have factorized form, namely, 

Qn,l{^l, Z2\Z3, Z^\... \z2l-\, Z2l\z2l+1 . . . Zn) = {Z\ - Z2){Z'i - Z4) . . . {Z2l-1 - Z21) W (Zj - Zk)- 

2l+l<j<k<n 

Qn,l{Zl, Z2\Z3, Z^l . . . \Z2i-l, Z2l\z2l+1 ■ ■ ■ Zn) (3.23) 

where the polynomials Qni{zi, Z2\z3, z^l . . . \z21-1, Z2i\z2i+i ■ ■ ■ Zn) depend only on differences of vari- 
ables Zj — Zk- Also they have the same symmetry properties as the function 

An,i{zi, -^21-^3, Zi^\ . . . \z21-1, Z2i\z2i+i ■ ■ ■ Zn) described above. The maximum powers for each variable 
may be readily calculated from the corresponding power of the polynomial Qni- 
Let us rewrite the properties 2 and 3 from the Section 2 in the following form 

Pn{Zi,. . . ,Zn-2,Zn-l,Zn-l±i) = (3.24) 
lim Pn{zi, . . . ; Zn—1: Zn ) -Pn-l{Zi, . . . , Zn-l) (3.25) 

As we shall see below together with our main ansatz (3.^0]) they completely fix the answer for 
Ph{z\., . . . , Z5) and P%{z\., . . . .,z^. Therefore we may expect that it also happens for any P„( Z\i ■ ■ ■ 1 Zn ) ■ 

Let us start our consideration from the analysis of the relation (3.^5]). First let us make a simple 
observation that it has an obvious solution corresponding to the case when all the variables Zj go to 
infinity § i.e. 

P„(oo,...,oo) = ^ (3.26) 

Now we can analyze corollaries of the relation (3.^) and our general ansatz (3.^Up. It can be 
seen that this is equivalent to the following chain of recurrent relations for the functions An^i 

^2fc+l,K^l5 Z2\zz, Z^\... \z2l-\, Z2l\z2l+l ■ ■ ■ Z2k, Oo)- 

^{^2fc+i,«+i(^i, ^2^3, 2:4! . . . |oo, Z2i+i\z2i+2, ^2fc) + pcrmutations} = 



1 

-A2k,l{^iy^2\z3,Z4\...\z2l+l,...,Z2k) foT l = 0,l,...,k (3.27) 



2 

2 

A-2k,l{^l^ ^2^3, Z4\.. . \Z21-1, Z2l\z2l+l ■ ■ ■ Z2k-1, Oo)- 

^{^2^,^+1(2:1, ^2^3,^41 • • • \oo,Z2i+i\z2i+2, ■ ■ ■ , Z2k-i) + permutatious} = 

^^2fc-i,«(^i)^2k3,^4| • • • k2Z+i, • • • ,^2fc-i) for / = 0, 1, . . . , - 1 (3.28) 



Since An^o are just constants below we will not write their arguments any more. 
The limits z„ — > cx) , . . . , zi — > cx) should be taken consequently 
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where permutations are taken only for the variables Z21+1, Z21+2, ■ ■ ■ Z2k for (3 .[271) and Z21+1, Z21+2, ■ ■ ■ Z2k-i 
for (3 .[281) in such a way that the variables within symmetric sets are not permuted. In the formula 
(3. |27|) it is also implied that 

A2k+i,k+i = 0. (3.29) 

Up to the moment it is not completely clear how to write down in a closed form all the equations 
for the functions An^k which follow from the relation (3.|2^. 

Let us write them down together with the equations (3.^) and (3.^8|) consequently for n = 
2,3,4,5,6. 

The case n = 2 

Since A; = 1 in this case we have only one equation in the set (3.|2^) 

A2,o - ^A2,i(oo,2;i) = ^^1,0 (3.30) 

which is obvious because as it can be seen from eqs. (3.|^3.|^), (3.7) and (3.|^) all the functions A 
here are just numbers, namely, 

1 

^1. = 2 

^2,0 = I, A2,l{zi,Z2) = ^ 

3 o 

The second property, namely, the formula (3.|2^) is also trivially satisfied because using eq. (3.|T^), 
(3.|T1) we have from (3.D 

P2{zi,zi + i) = A2fi + A2,i{zi,zi + i)G{-i) = 

A2,o - 2A2,M,z, + i) = (3.31) 
So the equation corresponding to the second property is just the last relation in formula (3.^). 

The case n = 3 

Since n is odd the first property (3.[25|) corresponds to the relations (3.[27|) with k = 1. Here we have 
two of them ^ 

As^l{Zi, Z2\00) = -A2,l{Zi,Z2) 

^3,0(2:1, 2:2 1 00) - ^1^3,1(00, 2;i|2;2) + A3^i{oo,Z2\zi)^ = ^^2,1(2:1,^2) (3.32) 
The second property (3.^) i.e. the equation 

^3(^1,^2,2:2 + 2) = (3.33) 
is equivalent to the following two relations 

^3,1(2:1, 2:212:2 + i) +7(^12)^3,1(2:1,2:2 + 212:2) = 

^3,0 + ^(2:12)^3,1(^1,^2 + 212:2) - 2^3,1(^2,2:2 + 212:1) = (3.34) 
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Substituting here ^3^1 defined in (3.7) and (3.||) {A^fi = 1/4) one can easily check that these two 
equations are satisfied. 



The case n = A 

As in the previous case the first property (3.^51) is equivalent to the set of two relations (see (3.^) 
for k = 2) 

-44,1(^1,^2^3,00) - 2^4,2(^1,^2100,2:3) = -A3^i{zi,Z2\z3) 
Aifi - ^I^A^^i^Zi, 00\Z2, Z3) + A4^^i{z2, 00\zi, Z3) + A4^^i{z3, 00\zi, Z2)^ = ^Asfl (3.35) 

The Property 2 (3.^4]) for n = 4 is equivalent to the following set of equations 

7{z23)A4,2{Zl,Z3\z2,Z3 + i) + 'y{Zi3)A^^2{Zl,Z3 + i\z2,Z3) = 

Ai^l{Zi, Z2\Z3, Z3 + i) - 2Ai^2{Zl,Z2\z3,Z3+i) = 
A4^l{zi,Z3\z2,Z3 + i) + -f{zi3)A4^i{zi,Z3 + i\z2,Z3) + « (2:23)^44,2 (^1 , 2:3 ^2 , 2^3 + = 

^4,0 + 0(213)^4,1(^1,^3 + ^1^2,23) + "(^23)^4,1(^2,^3 + ^1^1,^3) - 2^4,1(2:3, 2:3 + z|2;i, 2:2) = (3.36) 

Again using eqs. (3.7) and (3.p!0|-3.pl]) one can check that the equations (3.|35|-3.|36D are satisfied. 

So we have established that, in fact, the formulae (3.7) with the definitions (3.^3.|TT]) provide the 
solution to the equations (3.|30D and (3.^) for n = 2, eqs. (3.^2]) and (3. [3^) for n = 3, eqs. (3.p5|) 
and (3.|36|) for n = 4. 

We may also go an opposite way, namely, demand all these equations, solve them in framework 
of the ansatz (3.^0|-3.^3D and then come to the result for P„. Let us mention that this way seems to 
be much easier than the direct calculation of the integral (3.|I]) as we did before. 

For n = 5, 6 we would like to try this new way of getting the result for EFP. 

The case n = 5 

First as above we write down consequences of the properties (3.^Sp and (3.^ip. The first one (3.p5|) 
is equivalent to the set of three relations which can be obtained from eqs. (3.|27|) by taking k = 2 

^5,2(2:1, 2^2 ^3, 2:4 1 OO) = ^A4_2( 2:1,^2 1 2:3, 2:4) 
A^^li^Zi, Z2\Z3, Zi, 00) - ^i^A5^2{zi, Z2\00, Z3\Z4) + ^45,2(2:1 , 2:2 | OO, 2:412:3) | = ^^4,1 (zi, ^2 ^3, 2:4) 

As^o - ^1^5,1(00,^11^2,2:3,^4) + v45,i(oo, 2:212:1, 2:3, 2:4) + 

+ ^5,1(00, 2:312:1, 2:2, 2:4) + ^5,1(00, 2:412:1, 2:2, 2:3)) = -^^4,0 (3.37) 



2 

The Property 2 (3.53) for n = 5 is equivalent to the following set of five equations 



^5,2(2:1,2:212:3,2:412:4 + i) + 7(2:34)^5,2(2:1,2:212:3,2:4 + i|z4) = 
7(2:14)^5,2(2:1,2:4 + ^122,2:412:3) + 7(^24)^5,2(2:1,2:41^2,2:4 + ^12:3) = 
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a{zu)A5,2{Zi,Z2\z3,Z4 + i\z4) - 2Ar,^2{Zi,Z2\z4 + i,Z4^\z3) + Ar,^i{zi, Z2\z3, Z4^, Z4 + i) = 



"(2:24)^5,2(^1, ^4^2, + «|2;3) + a{z34)A5^2{Zl, Z^lz^ + i, Z3\Z2) + A5^i{Zi,Z4\z2,Z3,Z4 + i) + 

7(^14)^5,1(^1,^4 + ^^2,^3,^4) = 

^5,0 + a{Zi4)A5^i{zi,Z4 + i\z2,Z3,Z4) + a{z24)A5^i{z2,Z4 + i\zi,Z3,Z4) + 

+ a{zu)A::,^i{z3,Z4 + i\zi,Z2,Z4) - 2A5^i{z4,Z4 + i\zi,Z2,Z3) = (3.38) 

The difference with above cases n = 2,3,4 is that now up to the moment we do not have a 
solution to these equations. Therefore we have to get it. Let us briefly describe how we do this. First 
we take the ansatz (3.pi|-3.[23|), namely, 

A,o = I (3.39) 



A r \ \ Q->,li.Zi,Z2\z3,Z4,Z^) to An\ 

A5^l{Zi, Z2\Z3, Z4, Zr,) = — (3.40) 

2:i3^142;i52;23 2^24 2^25 

^5,2(^1,^21^3,2^4^5) = ^5,2(2^1,2:212:3,2:412:5) ^g^^^^ 

2:132:142:152:232:242:252:352:45 

and demand that the polynomials i(-2i, -2212:3, ,24, ,25) and 2(2:1, 2:2! ^3, -2412:5) have the same sym- 
metry properties as the functions 1 and ^5^2 and the same maximum powers of the variables 
Zi, . . . ,Z5 as their denominators in the r.h.s. of (3.^0|) and (3.|41|) respectively. Namely, we write 

Q^^^{zi,Z2\z3,Z4,Z5) = C5^i{ii,i2\i3,H,H) z\^z^z^2^lzf (3.42) 

< ii, ^2 < 3, < Z3, 24, ^5 < 2 
ii + . . . + Z5 < 6 

Q5^2(2:i, 2:212:3, 2:41^5) = C5^2{ii,i2\i3,i4\H) ztzi^zi'zfzi" (3.43) 

< ^1,^2,^3,^4 < 3, < ^5 < 4 
Zl + ... +^5 < 8 

with the coefficients C5^i(ii, 22^3, "^4, "^5) and C5_2(ii, ^2^3, "^4^5) which have evident symmetry proper- 
ties 

C5,i{ii,i2\i3,H,i5) = C5^i{i2,ii\i3, 14,15) = ^5,1(^1, ^2|o-(^3), o-(«4), (^(^5)) (3.44) 
where a is any element of the permutation group of three elements 5*3 while 

05,2(11, i2\i3,H\i5) = ^5,2(^2, ^1^3, ^4^5) = C'5,2(^l,^2K4,^3K5) = C'5,2 (^3, ^4^1 , ^2 ^5) (3.45) 

In accordance with eq. (3.|19D we also demand the translational invariance of these polynomials 

Q^^-i^{zi + a,Z2 + a\z3 + a,Z4 + a,Z5 + a) = Q^^-^^^zi, Z2\z3, Z4, z^) (3.46) 
Q.2izi + a,Z2 + a\z3 + a,Z4 + a\z5 + a) = Q^o{z\,Z2\z3,Z4\z5) (3.47) 
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These conditions allow us to fix some of the coefficients Cs^i and 6*5,2 in (3.^2|-3.p3D. 

The rest of these coefficients may be fixed by satisfying the equations (3.|38|). One should say 
that in spite of the fact that it turns out to be an over determined system it has a solution and this 
solution is unique. Also it is interesting to note that after the equations (3.|38D are satisfied all the 
residual equations (3.[37|) are satisfied automatically (!). 

Let us show the result 

^5,2(^1? ^2|^3, ^412^5) = 

^(^15^25 - ^){Z35Z45 - 1) (54,2(2^1, Z2\z3, Z^) + ^(^12 + ^){zli + 4)(5 - 2:132:24 - 2:142:23) (3.48) 

^5,1(^1; ^2|^3, ^4, ^5) = 

2 (^152:25 — 1) <54,2(^l5 ^21^:3, 2:4) + -(2:14^24 — 1) ^4,2(^15 ^2|2;3, 2:5) + -{Z13Z23 — I) Q4^2{zi, Z2\Z4, Z5) 

- ^^e^l^^ + ' 2)(^i5^25 + + (4 - 2)(2:i42:24 + + {zl - 2)(2:i32:23 + ^)} (3.49) 

Now using eqs. (3.^3.|4l|) and substituting these formulae into the ansatz (3.|20|) we get the 
answer for P5 in the inhomogeneous case. 

The case n = 6 

We proceed in the same line as for the case n = 5. The Property 1 i.e. eq.(3.^5D or equivalently the 
set of the recurrent relations (3.^8]) at A; = 3 gives us three equations 

-46,2(^1, 221^3,24125, 00) - ^yl6,3(2i,22|23, 2412:5, 00) = ^ ^5,2 (21, 22 I23, 24 12:5) 



2 ^' ' 2 

2 



^6,1(21, 2212:3, 2:4, 2:5, 00) - ■^<{ ^6,2 (21, 2:2 1 00, 2:3 1 2:4, 25) + ^6,2(21,2:2100,2:4123,25) + 



(21, 22I0O, 2:512:3, Z4)| — -v45,i (2:1, 2:212:3, 2:4, ^5) 
^6,0 - ^1^6,1(00,2:112:2,2:3,2:4,2:5) + ^6,1(00,^21^1, 23, 24, 25) + v46,i(00, 2:312:1, 2:2,^4, 25) + 

+ ^6,1(00, 2:412:1, 22,23, 25) + ^6,1(00, 2:512:1, 2:2, 2:3,^4)! = ^^50 (3.50) 
The Property 2 (3.^4|) for n = 6 produces seven equations 

A,2(2i, 2212:3, 2:412:5, 2:5 + i) - 2^6,3(21,2:212:3,2:412:5,2:5 + i) = 

7(245)^6,3(21,22123,251^4,25 + + 7(2:35)^6,3(21, Z2|23, 25 + Z|2:4, ^5) = 
^6,2(21,22123,2:512:4,25 + Z) + 7(235)^,2(21, 2212:3, 25 + Z|Z4, 25) + 

+ ^(2:45)^6,3(21, 2:212:3, 2:512:4, 2:5 + i) = 
0(2:35)^6,2(21, 2:212:3, 2:5 + i|2:4, 25) + 0(245)^6,2(21, 2:212:4, 2:5 + z|2:3, 2:5) - 
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-2AQ^2{Zi,Z2\z5,Z5+i\z3,Z4,)+A(i^l{Zi,Z2\z3,Z4,Z5,Z5 + i) = 



7(2^25)^,2(21,^51^2,^5 + ^1^3,^4) +7(215)^,2(21,25 + ^1^2,25123,24) = 

^6,1(21, 251^2, 23, 24, Z^+i) + 7(2:15)^6,1(^1, Z5 + i\z2, Z3, Z4, 25) + 

"(225)^,2(21, Z5\z2, zr, + i\zs, Zi) + 0(2:35)^6,2(21, 2:5123, 2:5 + i|2;2, 2:4) + 
+ 0(245)^6,2(21,2:512:4,2:5 + ^1^2,23) = 

^6,0 - 2^6,1(2:5,2:5 + i|zi, 2:2, 2:3, 2:4) + 0(2:15)^6,1(2:1,2:5 + ^1^2,23,24,2:5) + 

«(225)-46,l(22, 2:5 + 2|2i, 2:3, 2:4, 2:5) + 0(2:35)^6,1 (2:3, 2:5 + 2|2;i, 2:2, 2:4, 2:5) + 

+ 0(^45)^,1(24,25 + ^12:1,2:2,2:3,2:4) = (3.51) 

In the same line as we did it above for the case n = 5 we can solve the equations (3.|50|) and (3.pl|) 
in framework of the ansatz (3.pT|-3.|2B|) 

v46,o = \ (3.52) 

A I \ \ _ <56,1 (21, 22 1 23, 2:4, 2:5, 2:5) /Q e;Q\ 

^6.ll2i, Z2|23, 24, 25, 2:6j — (^O.OOJ 

2132142152:162:232:242:252:26 

^,2(21, 2212:3, 2:412:5, 2:6) = ^6,2(21,22123,2:412:5,2:6) {^M) 

2132142152:162:232:242:252:262352:362:452:46 
A / \ \ \ Q 6,3(^1, Z2\Z3, Z4\Z5, Ze) ro 

^6,3(21,221^3,24125,2:6) = (3.55) 

2l3 2l4 2i5 2:16 2:23 2:24 225 226 235 236 245 246 

As above we demand that the polynomials (^6 1(21, 2:212:3, 2:4, 2:5, 2:6), (562(21,22123,2:412:5,2:6) and 
(^6,3(21, 2:212:3, 2:41^5, 2:6) have the same symmetry properties as the functions Aq^i, Aq^2 , ^6,3 respec- 
tively. Analogous to eqs. (3.^2|-3.^) we write 

^6,1(^15 ^2|23, 24, 2:5, 2:6) = ^ (^6,1(^1,^2^3,^4,^5,%) 21^2:2^ 2:3^ 2:4" Zg^Zg*^ 

< ii, 22 < 4, < 13, ii, 15, iQ < 2 
ii + ... + iQ<8 

_ (3.56) 

^6,2(21, 22123,2:412:5, 2:6) = J2 Cfi,2iiui2\i3,i4\i5,i6) z[^4^4'4^4'4'' (3-57) 

< ii,i2, • • • ,^6 < 4 
ii + . . . + i6 < 12 

Q6,3(2i, 22 123, 2:412:5, ze) = Ce,3{h,i2\i3,iA\i5,H) 4'44'4*4'4^ (3.58) 

< ii,i2,...,i6 < 4 
n + . . . + ^6 < 12 

with the coefficients C6,i(ii, ^2^3, ^4, ^5, ^e), (^6,2(^1, 'j2K3, ^4^5, ^6) and C6,3(ii, ■j2K3, ^4^5, ^6) which sat- 
isfy the symmetry conditions 

^6,1(^^1,^2^3,^4,^5,^6) = C6^i{i2,il\i3,H, 15,16) = C6,l(n,^2|cr(^3),Cr(^4),0-(^5),0-(^6)) (3.59) 
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where a is any element of the permutation group of four elements 6*4 

C6,2(H, ^2^3, ^4^5,^6) = C'6,2(ii2, ^1 1^3, U^S, ^^e) = C'6,2(^l, ^2^4, ^3^5, ^e) = C'6,2(^l, ^2^3, ^4^6, ^s) = 

= C6,2{i3,U\ii,i2\i5,i6) (3.60) 

C6,3(^1,«2K3, ^4^5,^6) = C'6,3(^2, ^1 1^3, ^4^5, = C'g.sI^I , ^2 1^4, «3 1^5, ^g) = C'6,3(«l , ^2 1^3, ^4^6, ^s) = 

= Ce^3{i3,u\ii,i2\i5,i6) = C'6,3(»ii, ^^2^5, ^6^3, u) (3.61) 
As in the previous case the translational invariance of these polynomials 

QQ^i{zi + a,Z2 + a\z3 + a,Z4^ + a,Z5 + a,ze + a) = Qe^i{zi, Z2\z3, Z4, Z5, zq) (3.62) 

Qe^2i^i + a,Z2 + a\z3 + a,Z4^ + a\z5 + a,Z(i + a) = Q^^^i^i, Z2\z3, Z4\z5, zq) (3.63) 
QQ^3{zi + a,Z2 + a\z3 + a,Z4^ + a\z5 + a,ZQ + a) = Qq^^{zi, Z2\z3, Z4\z5, zq) (3.64) 

allow us to fix a lot of the coefficients Cg^i, Cq^2 and Cg^s. 

Then with the help of computer we solve the over determined system of the equations (3. |5l|) . It 
is left to check that all other equations (3.|50D are satisfied automatically as in the case n = 5. 

The result looks as follows 

Qe,^^^!^ ^2\Z3, Z4\Z5, Zq) = 
6^ (54.2(^1, ^2^3, ^4)Q4,2(^1> ^2^5, 2;6)Q4,2(^3, ^4^5, Zq) + {zI^+A){zI^+A){zIq+A) A6,3(^1, ^2^3, ^4^5, ^e) 

(3.65) 

Q6,2i^l^ Z2\Z3, Z4\Z5, Zq) = 
6^ — — — 2 2 

— Q4,2(^l5 ^2^35 Z4)Q4^l{Zi, Z2\Z^, Zi^Qj^^^{z3, Z/^Z^, Z^ + {z^2 + ^{"^3/^ + 4) Ag, 2(^:1, Z2\Z3^ Z^lz^, Zq) 

_ (3.66) 

Q 6,l{Zl, Z2\Z3, Z4,, Z^, Zq) = 
^ (^13^23 - 1)(^14224 - l)(2l5^25 " 1)(^16^26 - 1) + (^12 + 4) A6,l(2;i, ^2|^3, Z4, Z5, Zq) (3.67) 



where 



A6,3(2;i, ^2|^3; Z^lz^, Zq) — g^gg (5 ~ ^13^24 " ^14^23)(5 — ^15^26 ~ ^16^25)(5 — ^35^46 ~ Z3QZ45) 

^ (4 - 26)(4 - 26)(4 - 26) + -^(4 + 4 + 4 - 16) (3.68) 



189000' ' 200 



^6,2{Zi, Z2\Z3, Z4\Z5, Zq) = | ^ A6,3 (^1 , ^2 1^3, ^4 1^5, ^e) + ^^^42 (^1 , 2^2 ^3 , 2:4) (4 + 1) } (^56 + ^3) + 

|y^(5 — 2132:24 ~ 2;i42;23)(4 " 3 + Zi5Z2e + Z1QZ25 + 2:35246 + 236245) — 

2^(4 + y)(%Uy)-^}(4 + l) (3.69) 
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Z2\Z3, 24, Z^, Zq) — 77— (2:13 2:14 2:23 2^24 2:56 + Z13Z15Z23Z25ZIQ + ... + Zi5ZiqZ25Z2qzI^) + 

oUU 

900 ^^^'^^^^^^^^^^ + 2:13 2:15 2:23 2^25 + ••• + Zi5ZieZ25Z26) — 
— 225!^^"'^^^^'^ ~'~ ^1*^^23)2;56 + (^132:25 + Zi^Z2z)z\q + . . . + (2:15^26 + 2:16 2:25) 2^34 1 — 
~ 300 1''^"^^^^^'^^"^''^^^'' ^ (2;i32;25 + ^152^23) + • • • + [zi^Z2& + ZifsZ2^)^ + ■^^zl2{zl4zlf^ + zl^zlt^ + zlf^zl^) - 

~ 2520 ^^^^^^^ ~'~ ^35^46 + ^36^45) + 1575 ^-"-^^^^^ ~'~ ~'~ ~^ ~'~ ~^ ^^^^ ~ 

- Y^(4 + 4 + 4 + 4 + 4 + 4) + ^^?2 + ^ (s.ro) 

By means of this solution, eqs. (3.|52|-3.|55D and the ansatz (3.|20|) for n = 6 we come (for the very 
first time (!)) to the resuh for Pg in the inhomogeneous case. 

Let us make a remark about some additional amusing factorization property of the polynomials 
Qj^i whose meaning is still to realize, namely 

Q3,i{zi,zi + 2i\z3) = ^{zi3 + if 

QA2{zi,zi + 2i\z3,Z4) = ^{zi3 + if{zM + if 

3d 

(54^i(zi,zi + 2z 1^3,^4) = ^{zi3 + iy{zu + if 

Q^2izi,Zi + 2i\z3,Z4\z5) = ^{Zi3+iY{zu + iY{zi5 + if Q^^{Z3,Z4\Z5) 

D 

Qr,^^{zi,zi + 2i\z3,Z4\z5) = ^{zi3 + if{zu + if{zi5 + if 
Qq3{zi, zi + 2i\z3, za\z^, zq) = ]:{zi3 + if{zu + ifizi^ + if{ziQ + if Q4 2(^3, -24^5, zq) 

D 

Qq2{zi, Zi + 2i\z3, Zi\z^, ^e) = 7(^13 + if{ziA + if{zi-> + if{zi(i + if 1(^3, ^4^5, ^^e) 

D 

QQ^i{zi,Zi + 2i\z3,Z4\z^,Z(i) = ^{zi3 + if{zu + if{zi5 + if{zi(i + if 
Looking at these formulae we can assume that 

1 

-n(^i.+o' (3.71) 

" J=3 







\z2l-lZ2l\z2l+l ■ ■ - Zn) 


Qn-2,l-l{z3, ZaI 




\Z2I- 


-lZ2l\z2l+l ■ ■ ■ Zn) 
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4. The structure of P{n) in homogeneous case 



As we have discussed in ||I0[ the formulae in the homogeneous case become more regular in terms of 
the alternating zeta values 

(a{s) = E = -Li.(-l) (4-1) 

where Lis(x) is the polylogarithm. The alternating zeta series is related to the Riemann zeta function 
as follows 

C(^) = Y^^Cais) (4.2) 

Unlike Riemann zeta function the alternating zeta series is regular at s = 1, it is Ca(l) = log 2. 
Let us mention that the answers for P(l), . . . , -P(4) can be obtained from the formulae (3.^3.^ by 
taking the homogeneous limit i.e. Zj and using the expansion of the function G (3.p^. 

The same can be also done for n = 5 and n = 6 hj means of the formula (3. |20|) , formulae (3. |2T|) , 
(3.p^), (3.|2B|) and the answers for the polynomials (3.p5|-3.^DD and (3.^3|-3.^^ for the cases n = 5 
and n = 6 respectively. 

For n < 5 we reproduce the known results, see |10| formula (1.16) in there. 

In the case n = 6 we discover a new result 

1 r 9244 22694 

P(6) = -|l-35Ca(l) + 322C„(3)-^Ca(5) + ^— Ca(r)- 2982Ca(9) 

- ^ ai) ■ c(3) + ^ ai) ■ c(5) - ^ usr - ^ ud ■ ar) 
+ ^ US) • as) + ^ ai) ■ CM ^ ^ US) ■ U7) - ^ u^)^ 
- ^ US) ■ U9) + ^ a5) ■ U7) + ^ a5) • U9) - ^ ur)^ 
. ^ ai) ■ US) • a5) - ^ a3)3 - ud ■ us) ■ U7) . ^ ■ u^)^ 
- as)^ • US) . ai) ■ US) ■ a9) - us)^ ■ ar) + us) ■ as)^ 

- ^ UD ■ C(5) ■ a9) + ^ ai) ■ C(r)^ + ^ US)^ ■ C(9) 

264082 ^ ^ ^ 188630 
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C(3)-a5)-C(7) + ^^Ca(5)^} (4.3) 



This support our hypothesis that all P{n) can be expression in terms of values of Riemann zeta 
functions at odd arguments, log 2 and rational coefficients. 
From this expression we can get numerical value 

P(6) = 7.068127533 ■ 10"^ (4.4) 

In [|T0| a numerical method was used for evaluation of P{Q). It is called Density Matrix Renormal- 
ization Group (DMRG). The results can be found in Table 1 in [|10]. In particular P{Q) = 7.05 ■ 10~^ 
with an uncertainty in the second digit after the decimal point. It is in a good agreement with our 
analytic result (4.3), (4.|D. 
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Looking at the formulae (1.16) of [|10] and the above expression (4.3) we can make a general 
conjecture for the dependence of P{n) on the alternating zeta series 



n-2 



Pin) = j-—^ E n [C(2 J + 1)]^^- (4.5) 

where all coefficients -B^" r„_2 rational and the sum is over non-negative integers ro, . . . ,rn-2 
which belong to the region U determined by the following two conditions 

n-2 
j=0 

E r, {2j + 1) < -^—^ (4.6) 

j=0 ^ 

Let us show how the non-zero coefficients -Bf"!., ^ , look like for the cases when we know the 
manifest analytic answer, namely, when n = 1, 2, ... 6 

5(1) = 1 
Sj^) = 1 = -1 

75(3) 1 0(3) /] d(3) o 

-Dq.O - J- -Dl,0 - -4 i^o,l - ^ 

r(4) _ 1 r(4) _ - ^'^^ r(^) - 

-°o,o,o — -°i,o,o — -LU -Do,l,0 — -°o,o,i — 

r(4) _ (4) _ 170 (4) _ 17 

-Dl,l,0 — ^ -°1,0,1 — -^0,2,0 — ^ 



r(5) _ 1 r(5) _ _<.n _ 281 (5) _ _1355 (5) _ 889 

-^0,0,0,0 ~ -'-'i.o.o.o ~ '^^ -°o,i,o,o ~ 2 -"0,0,1,0 ~ g -"0,0,0,1 ~ g 



r(5) _ ion r(5) - r(5) - r(5) - ^''^ r(5) - fi7Q 

-"1,1,0,0 ~ -"-"'J -"1,0,1,0 ~ ^ -"1,0,0,1 ~ ^ -"0,1,1,0 ~ -"0,1,0,1 — u'y 

r(5) _ qofi r(5) _ 970 
-"0,2,0,0 — — -J^O -Do,0,2,0 — ^ 

0(6) _ . 

-"0,0,0,0,0 — 

r(6) _ qp- r(6) _ 099 r(6) _ r(6) _ 22694 (6) _ 

-"1,0,0,0,0 ~ 'J"-' -"0,1,0,0,0 ~ -J^^ -"0,0,1,0,0 ~ g -"0,0,0,1,0 ~ g -"0,0,0,0,1 ~ ^yoz 

(6) _ 3920 (6) _ 369908 (g) _ 263816 (g) _ 323344 

-"1,1,0,0,0 — ^ -"1,0,1,0,0 — -"1,0,0,1,0 — ^ -"1,0,0,0,1 — ^ 

.6) _ 3458 (6) _ 933702 (e) _ 2627842 (e) _ 235963 

-"0,1,1,0,0 — -"0,1,0,1,0 — ^ -"0,1,0,0,1 — -"0,0,1,1,0 — ^ 

(6) _ 368564 (6) _ 28784 (e) _ 751592 (e) _ 644987 

-"0,0,1,0,1 — ^ -"0,2,0,0,0 — ^ -"0,0,2,0,0 — ^ -"0,0,0,2,0 — ^ 

(6) _ 538496 (6) _ 1143268 (g) _ 1737148 (g) _ 528164 

-"1,1,1,0,0 — -"1,1,0,1,0 — ^ -"1,1,0,0,1 — -"1,0,1,0,1 — ^ 
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(6) _ 264082 (6) _ 653296 (g) _ 924287 (g) _ 124082 

-DQ, 1,1, 1,0 — g -Ol, 0,2,0,0 — g -^1,0,0,2,0 — g -^0,1,2,0,0 — g 

(6) _ 163324 (6) _ 1737148 (g) _ 264082 

-D0,2,l,0,0 - -^0,2,0,1,0 - ^0,2,0,0,1 " ^ 



-^0,3,0,0,0 — -^0,0,3,0,0 — ^ (4-7) 

For us it was a bit surprising that two of coefficients B^^^ appeared to be zero, namely, 

B[%i,i,o = 5iS,o,o,o = (4.8) 

It means that the structures Ca(l) ■ Ca(5) ■ Ca(7) and Ca(l) ■ Ca(3)^ do not appear in the final answer for 
P(6). Meanwhile the term Ca(3)^ survived with the non-zero coefficient -Bqsooo which we expected 
to be zero. 

The most evident conjecture that we can make looking at the formulae (4.7) is as follows 

4:ol..,o = 1, (4.9) 

f n \ 

where is the binomial coefficient. 

y m J 

The next our conjecture is less trivial. As appeared the coefficients B might satisfy some equa- 
tions. One of them has the following form 

E = (-i)'f"7) (4-11) 

ro+ri + ...+rn-2 =P \ ^ / 

where p is some fixed positive integer. The expression (4.|TTp can be easily verified for the first 
coefficients given by the formulae (4.7). For example, when p = 1 we get the following equation 

<ol..,o + <il..,o + • • • + <ol.,i = -n + 1 (4.12) 

where we have already made the conjecture (4.|TDp for the first term. We believe that there should 
be more equations like (4.[Tl|) which probably provide the rigorous expression for the coefficients B. 

We have pointed out in our previous work ||T^ that since ^0(1), Ca(3), Ca(5), . . . are very likely 
different irrational (or even transcendental) numbers i P{n) seem to be different irrational (or tran- 
scendental) numbers as well. This means that P{n) does not satisfy polynomial recursion relation 
with respect to the distance. 

In order to clarify the structure of the formula (4.^ let us formally replace all Ca(2j + 1) by one 
complex variable x. This will define a new function P{n,x). 

Then we can calculate P{n,x) using our conjectures (4.|^), (4.|lTD and some properties of the 
binomial coefficients. Namely, 

^(■^■^)-g(-^r )^ ,4.13) 



'^It was proven by Apery ||T^ that ^(3) is irrational. Then Rivoal |jT^ proved that one of the nine numbers 
C(5), . . . , C(21) is irrational. One of the most recent theorem proved by Zudilin says that one of the four values 
C(5), . . . , C(ll) is irrational. (See also the paper by D. Zagier and the paper by Yu. Nesterenko |pl|). 
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where 



1 ± ^/l - Ax ^ 

A± = (4.14) 



In particular when x ^ 1 one gets 

P{n,l) = { 0, n = 3fc + 2 (4.15) 



1, 


n = 


6k or 


0, 


n = 


3k + 2 


-1, 


n = 


■6k + 3 



Note that the alternating zeta values approaches 1 as the argument gos to infinity 

Jim Ca(s) = 1 (4.16) 
Another nice form of the formula (4.p!3|) can be obtained by substitution 

^ = . ^2 (4.17) 
4 cosh a 

Then 

Pin, \ ) = sinh[(n + l)a] ^^^^^^ 
4 cosh a 2"(?7, + 1) cosh asinha 



We see that when a tends to zero 

^) = i 
a^o ^ ' 4 cosh^ a 2 



hmPin, ^ _ ) = - (4.19) 



Let us remark that this result appeared to coincide with the limiting formula (3.^). We do not 
know if this is accidental or there is a reason for this. 

Let us briefiy discuss the generating function for the values P{n) 

oo 

^(y) = (4-20) 

n=0 

where -P(O) = 1 by definition. Taking into account the conjectures (4.^) and (4.p!0|) we can easily get 
the two first terms for the generating function "^{y) 

*(,) = -i!l(i^ + ^^,„2 + ... (4.21) 

y Hi-yr 



As we discussed in and [jl0| we expect that for n ^ 1 

P{n) ~ e"'^"' (4.22) 
If we substitute it formally into eq. (4.[20|) with y = then we can expect that 

oo 

^(e") ~ ^(m) = ^6"'^"'+"" (4.23) 

n=0 

SO that the function \E'(m) satisfies the functional equation 

— ~ 77/ 

^(u) + ^i-u) - 1 = es{-,e-^) (4.24) 
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where 6^ is the third Jacobi theta function with the nome e"**. So we can expect that the generating 
function (4.|20|) may be related with the eUiptic functions and have some automorphic properties. 

Another possibihty is to put into the r.h.s. of eq. (4.[20|) the values P{n, x) given by (4.|l3|) instead 
of Pin) then 



^,{x,y) = Yy'^Pin.x) = / In ^ V + V^ ^ 4a: 



Another generating function looks much simpler 

oo 

^\x,y) = ^(n + l)i/"P(n,a;) 



n=0 



1 — y + xy'^ 



(4.25) 



(4.26) 



5. Discussion 



The main point of this communication was to consider the emptiness formation probability in the 
inhomogeneous case. The basic advantage of the inhomogeneous case is that we have more parameters 
at our disposal and relation to the qKZ. We derived three general properties of Pn from qKZ, which 
appeared to be extremely useful. From our experience of the direct calculation of P„ with n < 4 we 
conjectured a general ansatz for P„ (3. |20|) . We established that the ansatz (3.|^) with the first and 
the second property from the Section 2 [ namely (2.^] completely fix the answer. The third property 
(2.|^) turns out to be a corollary. This observation allowed us to evaluate P5 and Pe very efficiently 
[ in the inhomogeneous case]. The homogeneous limit of P5 reproduced the expression, which we 
obtained in paper by going through hard and long computations. That time we were not sure 
that we will be able to calculate Pq at all. Now it is possible to do it very quickly just by taking the 
homogeneous limit of our result in the inhomogeneous case. In the next publications we are planning 
to prove our general ansatz (3.pOD and evaluate P{n) for arbitrary n. Actually we think that qKZ 
approach is so powerful that we will be able to evaluate any correlation function in the Heisenberg 
XXX model and to show that it has a structure similar to P(n) 

[f] I 

F{zi, Zn)eu...,e„ = Xli^n,!""'" IT G{z2j-i - Z2j) + permutatlous} (5.1) 

1=0 ' j=l 



where ej = ±1 

P(Zl, . . . , Zn)ei,...,e. 

and 




We may also expect that similar to the functions 
rational of their arguments zi, . . . ,Zn also. 



n 

= (GSin^'^IGS), (5.2) 

l±al 
2 

A„,^i = A+'f ■+ from (3.|20|) the functions A'^y'" are 
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